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Abstract 

An existence theory is established for a coupled non-linear elliptic system, known as “vortex 
equations”, describing the fractional quantum Hall effect in 2-dimensional double-layered electron 
systems. Via variational methods, we prove the existence and uniqueness of multiple vortices over 
a doubly periodic domain and the full plane. In the doubly periodic situation, explicit sufficient 
and necessary conditions are obtained that relate the size of the domain and the vortex numbers. 
For the full plane case, existence is established for all finite-energy solutions and exponential decay 
estimates are proved. Quantization phenomena of the magnetic flux are found in both cases. 
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I. INTRODUCTION 


In recent years, there has been an increasing interest in the study of held theoretic models 


governed by Chern-Simons theory 


17|. These models generally lead to very complicated 


systems of non-linear equations, which have presented new challenges to mathematical an¬ 
alysts. Particularly, in condensed matter physics, Chern-Simons theory can be used to 
describe the fractional quantum Hall effect (FQHE) 18l-l33|. Ichinose and Sekiguchi 251 


established an effective theory for topological excitations in the generic (m, m, n) Halperin 
state in Chern-Simons gauge held theory, which describes the FQHE in double-layer 2- 
dimensional electron systems. Under the radially symmetric ansatz, Ichinose and Sekiguchi 
25[, obtained a system of Toda-type equations, limited to single soliton conhgurations, and 
studied the numerical and qualitative behavior of soliton solutions. 

In section 2, from the Lagrangian of the model, we get a coupled gauged Schrodinger 
equations governing the bosonized electron helds in the upper and lower layers, as well as, 
a set of constraint equations arising from the electromagnetic vector potential. In the static 
case and by ignoring long-range inter-layer and intra-layer Coulomb repulsions, we reduce 
the afore-mentioned equations to a hrst-order system via a “hrst integral”. The reduced 
system is of the BPS type (named after the seminal works of Bogomol’ny 3^ and Prasad- 


Sommerheld 


35|). We give an energy lower bound in terms of the total hux and show that 


all hnite-energy solutions over the full plane are topological in our context. The BPS system 
is then transformed into a coupled non-linear elliptic system or “vortex equations”, which is 
a general form of the Toda-type equations of Ichinose and Sekiguchi 25| . We present three 
sharp theorems establishing the existence, uniqueness, and exponential decay of the vortex 
solutions, in addition to, quantization of the magnetic flux. 

In section 3, we consider vortex solutions over a doubly periodic domain. Via a varia¬ 
tional principle and a direct minimization problem, we use a weak compactness argument 
to prove the existence of vortices. We establish necessary and sufficient conditions, as well 
as, the uniqueness of the solutions. By a simple integration, we get quantized magnetic 
flux formulas. The necessary and sufficient conditions for existence give limits on the vortex 
numbers in terms of the domain size. 

In section 4 and 5, over the full plane, we establish the existence, uniqueness, and expo¬ 
nential decay estimates of all finite-energy solutions. Using a Choleski transformation, we 
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find a variational structure of the elliptic system. Through the weakly lower semi-continuity, 
Gateaux differentiability, and strict convexity of an energy functional, we prove the existence 
and uniqueness of vortex solutions. Additionally, we exploit the decay estimates to exhibit 
quantized magnetic flux formulas. 


II. BPS AND VORTEX EQUATIONS 


The Lagrangian describing the FQHE in double-layer 2-dimensional electron systems 251 
is composed of two terms, the matter term and the Chern-Simons term Lcs-, 




^cs, 


( 1 ) 


where 


= i'lp^ido - ittQ - ittQ )V’t + i'lPiido - ia^ + iaQ )'ipi (2) 

^cs = + ^cs{%) 

= + ^a-d,a ^. ( 3 ) 

There are two bosonized electron helds, which represent electrons in the upper and lower 
layers, denoted by V’o- (o' = 1,2 or respectively. M is the mass of the electrons and 
p and q are nonzero real numbers, and a~ are scalar potential helds corresponding to 
the U{1) ® U{1) local symmetry, p, z/, and A take on the values 0,1,2 and is the 
antisymmetric metric tensor. Denote ejk = e^jk- The gauge-covariant derivatives are dehned 
as 

dJ'*' = dj — iOj =F ia~ -\- ieAj, j = 1,2, (4) 

where the external magnetic held, B, is directed along the 2 ;-axis and in the symmetric 

B 

Coulomb gauge A,- =- ejkXk- We denote the temporal derivatives by 

j 2 

Dq = Oq iog T • (^) 

The term V{'ijj^,'ijji) represents the interaction between electrons like the Coulomb repulsion 
and short-range four-body interaction. is given by 

V(i’t, i’-i) = ^ (^tV't + Ml) ^ (Mf - Ml) ^ ^> (6) 
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and more specifically represents the potential between bosonized electrons, l/’t) 

the long-range inter-layer and intra-layer Conlomb repnlsions. The corresponding Enler- 
Lagrange eqnations of the Lagrangian density ([T]) are 


iDlip^ = 


{D\dI + DlDl)'^;^ 


2M 


it nt^ 


_ 2^ — — 


+ 




4vreo j 
1 


+ ^ 


\x — x'\ 

{DiDi + 


\X — X‘ 


2j9 _ _ — — 

+ \x — x'p 


+ , 
dvreo 

J \ k 

Fi2 = 2p (V'lV’t + Mi) 
Fi 2 = 2q {Mi - 

- ‘m 


ho = i— 
M 


17 - ■ P 

(i^i^lPi - 

= ‘m 



X — X' 


'dx' 


(7a) 


(7b) 


(7c) 

(7d) 

(7e) 

(7f) 

(7g) 

(7h) 


where and — 9iya“. We recognize eqnations fITal) and fiTbl) as 

a conpled ganged Schrodinger eqnations governing the bosonized electron helds in the npper 
and lower layers, respectively. Eqnations flTcj) and flTdl) are the Chern-Simons constraints. 
As a resnlt of not treating the electromagnetic potentials a'^ = (a+) and a~ = (a“) as 
background helds, we get the constraint equations 

In this paper, we will assume that the long-range inter-layer and intra-layer Coulomb 
repulsions are negligible and focus on the static case of the Euler-Lagrange equations. With 
these assumptions, and letting F^jy = d^bu — di^bfj_ — eB and F^^y = d^bi, — di,b^ — eB represent 
the rate of change of the magnetic helds, where b^ = + a~ and b^ = a'j^ — a~, the system 
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fl7a|) - fl7h|) is equivalent to 


+ ^{p+ ?)iV'tiVt+ ^{p - 
Mt = + ^(p - 

Bi2 = 2{p + q)\'il)^\‘^ + 2{p - q)\'il)i\‘^ - eB 
Bi2 = 2{p - q)\'ip^\‘^ + 2{p + q)\'ijji\‘^ - eB 
dibo = dibo = -^{P - (^aq) - iJaD^i’a) 

d2bo = ^(p + ('0(7^1 V’fT - 'IpaD^i’a) 

d2bo = ^(p- - 'tp^Dl'tp^) . 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 

(8f) 

(8g) 


Above, we use the summation convention on the indexes j and a and dehne e^i- = 1, = — 1. 

The system (I5c))-(|8gp is similar to the non-relativistic Chern-Simons theory due to Jackiw 


and Pi 


26l | and discussed by Yang [l7j . Recall that for any complex-valued functions V’ and 
r], we have the identity 




Using this identity we get 


djbo = 


M 


{p - ejkq)dj\il)^\^ - (u'JV’t - 
,t 


( 9 ) 


( 10 ) 




1 

M 


-Xp + ^jkq)dj\Xi? - 

+ {D]Xi + i^jkOA^ 

and a similar equation for djbQ. Introduce the operators as follows, 

D^. = dI ± iDI and = D\ ± iD^. 


Using flTT|) in fl8a|l and flSbjl we get 


boX^ = -^ (^+^- + Bi 2 - eB^ 'ip^ 


boXi = 


2 2 
— {p + q)\XXA^ + 

'll nl 


2M 


DXDi+B^2-eB)xP^ 


( 11 ) 


(12a) 


(12b) 


2 2 

-^{p-q)\XXAi +j^{p + q)\Xi?Xi- 
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As a consequence of identities (fTOj) to fll2bD . we obtain a “first integral” of the system 
fl8B-([8^ of the BPS type, 


-OlV't ~ (-^1 ~ = 0 (13a) 

= {D\ - = 0 (13b) 

Bi 2 = 2{p + g)|V’tP + 2(p - g)|V'^P - eB (13c) 

Bu = 2{p - g)|V’tr + 2(p + g)|V'tr - eB (13d) 

bo = ^(p + g)lV’tP + + ^ (13e) 

^ + + + (13f) 

Consider the total energy [25[, 


^ = 5F / I ^ " *-DaV>.,r" + eB (IV., 1 ^ + IViP) - eB (|V.,,„|^ + IVi.oP) [*, ( 14 ) 


where = VpT^ and = ■\/p/2 are the held conhgurations for the ground state of 
the fractional quantum Hall effect and p is the average electron density. Additionally, the 
ground state for the fractional quantum Hall effect requires the hlling factor u to satisfy 


2Trp TV 

eB p 


(15) 


As in 


25l | and using equations fll3cp and fll3dp . we dehne the total Chern-Simons hux by 


^cs = 


1 

2 




dx. 


(16) 


Hence, we may write the following energy lower bound 


r. eB ^ 

- 4Mp 


(17) 


and note that solutions to the “self-dual” equations (I13ap and (113bp correspond to lowest 
energy conhgurations. 

We are interested in hnite-energy solutions of the system fll3ap - (ll3fj) . Using fll3ap and 
fll3bD . the hnite-energy condition reduces to 


E = 


eB 


+ - (l^toT + I'^poH) dx < oo. 


(18) 
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In particular, over the full plane and to guarantee that the amplitude of the electron 
helds are equal to that of the ground state, we see that the hnite-energy condition imposes 
the boundary condition: 

^ iV'poP = ^ and -)■ iV'por = | as |a;| -)■ cx). (19) 

Solutions satisfying flT^ are called topological. Therefore, in our context, all hnite-energy 
solutions over the full plane are topological. 

To arrive at the coupled non-linear elliptic system of interest (i.e, “vortex equations”) 
we hrst introduce the complexihed variables d = di+ 182 , (3 = bi + ib2, and A = Ai+ iA2. 
Away from the zeros of the held from the self-dual equation fll3ap . we obtain the system 


ddlnlpj^l = id (3 — iedA 

(20a) 

dd\n \ip^\ = —idl3 -1- iedA. 

(20b) 

Summine fl2UaD and fl20b). and using; A = dd = dd. we get 


A In -0^p = i{df3 — d(3) — ie{dA — dA). 

(21) 

Let us note that 


8(3 = div(/5) — i{dib 2 — d 2 bi) and dA = —iB. 

(22) 

Hence, 


AlnIV’tP = + q)\ip^\^ + A{p - q)\ip^\^ - App. 

(23) 

In a similar manner, we mav use the self-dual eauation fll3b) to obtain 


AlnlV'tT = 4(p-g)|V'tr + 4(p + g)|^/'^r -4pp. 

(24) 

Let us denote the set of zeros of the helds and tp]^ by = {pi,. 

. ,PArJ and = 

{gi,..., ^^ 2 }, respectively. Note that the zeros of the helds and 

integer multiplicities. Dehne the functions Ui and U 2 by 

are discrete and of 

Ml = In ^ — In p and M 2 = In t/’j, ^ — In p . 

(25) 


Then, using the transformation x 1 —)■ \/ppx, we arrive at the vortex equations, 

Ni 


Ami = 4 (fciie“^ + - 1) + An ^ 6p. (x) 

4 — 1 

(26a) 

J — A 

N 2 


Am2 = 4 (fc2ie“^ + ^ 226 “^ - 1) + An'^dgj(x), 

(26b) 

t=i 
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defined for all x G where K = (kij) is the matrix 


K 


1 

P 


( p + q p-q 
p-q p + q 


( 27 ) 


and 6 p{x) is the Dirac distribution concentrated at the point P. The positive integers Ni 
and N 2 are called the vortex numbers. In this paper, we are only interested in the coupled 
system, for which p not equal to q. 

We now state our main existence, uniqueness, asymptotic behaviour, and quantized mag¬ 
netic flux theorems. In what follows, the notation |r2| and \K\, denotes the size of the domain 
and the determinant of the matrix K, respectively. 


Theorem II. 1 Consider the coupled non-linear elliptic system fl26ap - fl26bj) . over a doubly 
periodic domain D. Let the matrix K, given by ([27]), be positive definite. Then a unique 
solution exists if and only if the condition 

{ 2ti 27 r 1 

\n\ > max {k 22 Ni - ki 2 N 2 ), {knN 2 - /c 2 iA^i) | (28) 

is satisfied. 


Theorem II.2 Consider the coupled non-linear elliptic system (I26al) - (l26bp over the full 
plane and satisfying the boundary condition 


ln(2) as |x| ^ 00 . (29) 

Let the matrix K, given by fl27l) . be positive definite. There exists a unique solution to the 
system, satisfying the boundary condition exponentially fast. More precisely, we have the 
following exponential decay estimate, 

(mi + In 2f + {u 2 + In 2f < (30) 

when |a:| is sufficiently large, e G (0,1) is arbitrary, C^ > 0 is a constant depending on e, 
and Ao = 4 min 12, y |. Additionally, 

+ iVwal" < (31) 

when |a:| is sufficiently large, 6 G (0,1) is arbitrary, and Cs > 0 is a constant depending on 











Theorem II.3 In both the doubly periodic and the full plane cases, considered in Theorem 
ED and Theorem \IL2. there hold the quantized magnetic flux integrals 


J Budx = —27ipNi and j Bi 2 dx = —27rpA^2, (32) 

where the integration is evaluated either over the doubly periodic domain Vl or the full plane 


III. SOLUTION OVER A DOUBLY PERIODIC DOMAIN 


In this section, we prove Theorem III. II We first establish necessary conditions for the 


existence of a solution and derive the quantized magnetic flux integrals of Theorem III. 31 
over a doubly periodic domain. Through a Choleski transformation, as in [l4, [^|, we hnd 
a variational principle. We establish a coercivity condition and via a direct minimization 
approach, prove the existence of a solution to the system fl26ap - fl26bj) over a doubly periodic 
domain. 

By a doubly periodic domain, hi, we mean a periodic lattice cell with a “gauge-periodic” 
boundary [l7|. In other words, our solutions are periodic over a cell domain modulo gauge 


transformations as introduced by’t Hooft 3^. In this concrete situation, we can identify hi 
with the 2-torus Q = 

There are functions Mq : hi —)■ M and Uq : U ^ R, whose existence and uniqueness (up to 
an additive constant) are guaranteed by Aubin in [37|, satisfying 


Amq = 


'W 


Ni 


dvr 5pj (x) and Auq 
i=i 


47rA2 

'W 


N2 


+ 4:71 ^ 6 , 




[X 


(33) 


i=i 


Let mi,M 2 : 12 —)• R be functions satisfying the system fl26ap - fl26bp . Define vi = ui — Mq and 
V 2 = U 2 — Uq on D. Then, 


Avi = 4(A;iie“o+’^i + ki2e^o+^^ - 1) + 
Av2 = 4(fc2ie“0+"l + ^ 226 “°+"^ - 1) + 


dvrVi 

W 

dvrW 

w 


(34a) 

(34b) 


Integrating fl34ap - fl34bp . we obtain the following linear system of equations in the unknowns 
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J^e'^o+'^^dx and 


ku / e“o+’'Mx + fci2 / e“o+’'2(ix = |fi| - ttAT^ 

Jn Jn 

k2i f e“o+’^Mx + fc22 f =\n\-7iN2. 

Jn Jn 

/ // 

Solving for e^o+^^dx and e^o+^^dx we obtain the necessary conditions 


|fi| 


TT 


e^o+-idx = V - (k22JVl - kl2JV2)^ =Vl>0 


1^1 


\K\ 


TT 


e^o+-2dx = L^- {knN2 - k2im)-f- = m > 0 


(35a) 

(35b) 

(36a) 

(36b) 


In 


f ‘2,71 27r 

(^22-^1 ~ ^ 12 -^ 2 ) 5 (^11-^2 ~ ^ 21 -^ 1 ) ^5 


(37) 


2 —-'\K\ 

In terms of the size of the domain, inequalities f|36ap and fl36bp give 
|r2| > max 

and the necessity condition of Theorem III. II is established. Note that the vortex numbers 
Ni and N 2 are constrained by the size of the domain. 

Integrating equations (I34ap and fl34bll over hi, and expressing the results in terms of iV'o-P, 
we get 

[ {2{p +q)\7p^\‘^+ 2{p-q)\7p^\‘^-2p)dx =-2 ttpNi (38a) 

Jn 

[ {2{p - q)\'ijj^\‘^ + 2{p + q)\^pl\‘^ - 2p)dx =-27ipN2. (38b) 

Jn 

Therefore, from fll3cp . fll3dp . and flT^ we arrive at the quantized magnetic flux formulas of 
Theorem III.31 


and 


Bi2dx = —2TTpN2- 


(39) 


/ Bi2dx = —2TTpNi 

Jn 

In its current form, the system fl34a|) - fl34bj) does not have a simple variational structure. 
However, when the matrix K is positive definite, a variational principle for the elliptic system 
can be found via a Choleski decomposition, i.e., there is a unique lower triangular matrix L 
such that K = LLJ. To this end, consider the Choleski transformation 


wi = and W 2 = 7^{kiiV2 - /c 2 l^'l)• 

V\K\ \B\ 


The system fl34al) - fl34bD becomes 


Atci = 


■q + \/\K\wi _|_ ^^12 ^Un+{\K\w2+k2i\^/\K\wi)/kii _ q 

VW\ 


Aw 2 = /^e^0+i\Kh2+k2l^/Whl)/l^ll _ (^2, 


(40) 

(41a) 

(41b) 
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where 

The corresponding functional I : H\ —)• M is 

I{w,,w-2) = ^ |l|Vw,p + l|Vw2p + (43) 

+ ^e«o+(|i^b2+fc2ixA^«>i)/fcii _ _ ^-2^2Ida;. 


Let us show that the above functional I satishes a coercive lower bound. We use the 
notation x IV^’^(M^), = IV^’^(M^), and || ■ ||i ,2 to denote the norm of H^. 

Decompose as the direct sum of M and the set H^, i.e. = M © H^. The set 
is defined as the collection of all w G such that J^wdx = 0. Let Wi = Wi + Wi and 
W 2 = W 2 + W 2 , where tci,ta 2 G M and wi,W 2 G From the necessary conditions fl36ap 

and fl36bp . we get 


\/\K\wi = ln(r 7 i) — In ^ y 

{\K\w2 + k2i^/\l^wi)/ku = HV2) “ _ 


Applying the above decomposition to the functional 


I{wuW2)-^^\\Vwi\\l = (^l-\n{T]i) + \n(^J ^ 


+ T ^ V 2 ( 1 - 111(772) + In 


p“o+(l-f^|«)2+fc2l\/pd4«l)/A:il^^ 


Note that kn > 0, since \K\ > 0. By Jensen’s inequality we obtain 


In 


In I / j > 7 ^ [ Modx + ln|D|, 

\Jn J l“l Jn 

X+(|ic|^2+fc2iv1^m)/fen^^^ f M"da; + ln|D|, 

/ |D| Jn 


and the coercive lower bound 


I{wi,W2) > ^X]||Vw;i||2-/3 


i=l 


(44a) 

(44b) 


(45a) 

(45b) 


( 46 ) 


is attained, where /9 is a positive constant independent of tci and W 2 . 
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From the coercive lower bound (H6|) . we conclude that the functional I is bounded from 
below. So it makes sense to consider the direct minimization problem: 


m = inf < I{wi, W 2 ) 


Wi,W2EH^j. (47) 

Let be a minimizing sequence of (1471) . i.e., choose functions wf and W 2 in H^, 

where fc = 1, 2, 3,... so that 

I{w’l, W 2 ) ^ m as k ^ 00 and I{wl,wl) > > ■ ■ ■ > m. 

From fl46l) and the decomposition wf = wf + wf where i = 1,2, we may write 

1 ^ 




2=1 

and conclude that {Vtcf} all belong in 
Recall the Poincare inequality 

[ fix)dx <C f \Vf{x)\^dx, f e [ f{x)dx = 0, 

Jfi Jq Jo. 

where C > 0 is a suitable constant. By Poincare’s inequality we obtain 

.. 2 2 

I(wW2) > xj; l|v™‘||l - fS > -/?, (48) 

^ i=l i=l 

for some suitable positive constants ads, z = 1,2, and conclude that the sequence {w^} all 
belong in L‘^{VL). Therefore, is bounded in H^. Moreover, and belong 

in H^, since 




< — taf II 2 —)■ 0 as k ^ 00. 


(49) 


After all, we are in a reflexive space and without loss of generality, we may suppose that 


{101,102) {w[°°\ 102°°'^) E H 2 weakly as k ^ 00. 


(50) 


Let us show that the sequences {wi} and {tcf} of real numbers are also bounded. We 
will need the Trudinger-Moser inequality of the form 


f{x)dx = 0 , 


(51) 


where Ci and C 2 are positive constants. 
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From fl44ap . we may write 


\/\I^w’[ = ln(r 7 i) - In 
Taking absolute value on both sides of (]] 

\^/\I^w^\ <|ln(r/i)| + 




in 


and enlarging we have 
1 


In 




( 52 ) 


(53) 


In 




Apply flSTj) to the last term in 


to obtain 


Iv^^il < \ Hvi)\ + 2 


In 


e<dx] +CiC 2 yA^\ 


Vw’ldx 


(54) 


Since the sequence {tCi} is bounded in H^, we conclude that the sequence {tCi} is bounded 
in M. In a similar way, we may conclude that {w^} is bounded in M. Without loss of 
generality, we may assume that 




as /c ^ oo, i = 1, 2. 


(55) 


Dehne for i = 1, 2. Therefore, from (l50D and fl55D we get that 

{w^,W 2 ) G Hi weakly as k ^ oo. (56) 

To conclude that is the sought out solution of the minimization problem 

flT7)) . and hence a solution to the system fl26al) - fl26bH . we appeal to the weak lower semi¬ 
continuity of I. The uniqueness of the solution follows directly from the strict convexity of 
the functional /, which can be shown by a direct calculation of the corresponding Hessian 
matrix. Therefore, Theorem III. II is proved. 


IV. SOLUTION OVER FULL PLANE 


In this section, via a variational principle, we prove the existence and uniqueness compo¬ 


nents of Theorem III. 21 Our approach follows the ideas by Jaffe and Taubes in 


m 


and Yang 


17l |. In contrast to the doubly periodic case, the vortex numbers are not constrained by 


the domain size. 

In order to establish a variational principle, we transform the system fl26al) - fl26bp . satis¬ 
fying the boundary conditions fl29H . to an equivalent system, which can then be view as the 
Euler-Lagrange equations of a respective functional. 
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Consider the background functions Mq, Uq : —)■ M, depending on a real parameter /i > 0, 

defined by 


Ni 


N2 


Uo[X) = 


— ln(l +/r|a; — pjl and 
i=i 


Uo[x) = 


+ fi\x - Qj 

t=i 


-2^ 


Note that Uq{x),Uq{x) < 0 for all a; G and 


Ni 


Auq — 4 




A^i 


^ (/i+|a:-p^f)2 ^ 


A^i 


Au'i = -iY, 


A* 


iV2 


^ (p+|a:-gj|2)2 ^ 


Let 


Afi 


9'o{x) 




Ar2 


and ^"(a;)=4^ 




^ (p+|x-pj|2)2 |a;_g^.|2)2- 

We note that Mq, Uq, g^, and gQ all belong to L2(M2). Dehne Ui = Mi — Uq and V 2 = U 2 — Uq 
on M2_ Hence, the system fl26al) - fl26bD becomes 


Am = 4(^fcne“o+"i + fci2e“o+"^ - + 17o 
Am = 4(^fc2ie“o+"i + k22e<^^^ - + g^ 


(57a) 

(57b) 


where ni(a:) —>■ — ln(2) and V 2 {x) — ln(2) as |x| —)■ oo. 

When K is positive dehnite, we again use the transformation fl40|) . Thus, the system 
fl57ap - (157bp is transformed to 


where 


An;i = _ i] 

VW\\ J 

I ^-^12 f u” + ^(\K\w2+K2i^/\¥\wi) _ 1 V - 

7 

An;2 = + ^2, 


+ 4) and ^2 — 2 + |^(iCii5(o - iC2i5'o)- 


(58a) 


(58b) 


(59) 


However, the functions hi and h 2 are not in L2(M2^ ^^d this is an undesired property for 
what follows. To correct this issue, we dehne hi = ui + ln(2) and h 2 = '^2 + ln(2). Using this 
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new definition, the system fl57ap - fl57bp becomes 


Avi = 

Ah2 = 4Qfc2ie“o+"i 


+ ^fci2e“o+"^ - 
+ ^^ 226 “°+"^ - 


+ 9o 



(60a) 

(60b) 


where Vi{x) —)■ 0 and ^ 2 ( 2 ^) —t 0 as |a;| -^ 00 . With the use of fl40|) . the system fl60ap - fl60bp 
is equivalent to 


Awi = ( 61 a) 

_|_ (^u”+(\K\w2+k2iy/\K\wi)/kii _ _j_ ^ 

v^V ) 

Aw2 = 2{ e<+(l^l«'2+fc2ixA^«^i)/fcii _ 1 ) + h2, (61b) 


where 

h. = and h 2 = 

Moreover, hi, ^2 G L^(M^), wi{x) —)■ 0, and W 2 {x) —)■ 0 as |a;| —)■ 00 . 

The system fl61ap - fl61bp is well-dehned and is the Euler-Lagrange equations of the func¬ 
tional I : H 2 —)■ M, given by 



{Kng'^-K2,g',). (62) 


J(wi,m;2) = 


+ 


+ i|Vw 2 p + ^e“o(e^"'^ - 1) \dx 
g“o (e(^21\/tfdl«l + |^b2)/fcll _ 1) 1 .(^ 2 : 

4 


1 
2 
2 kii 

W\ 

hi 


(63) 


VW\ 


Wi + (^2 - ‘ 2 )w 2 >dx. 


To prove the existence of a solution of the system fl61ap - fl61bp . it is sufficient to show 
that the functional I, dehned by (ES]), attains a unique interior critical point in some open 
ball in Hi. It is straightforward to show that the functional I is strictly convex. Gateaux 
differentiable, and hence weakly lower semi-continuous. To prove the existence of an interior 
critical point, we just need to show that the Gateaux derivative of / satishes a coercive lower 
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bound. We begin by rewriting the functional I in the form 


+ - 1 - ^ik2l^A^\wl + 

|ff| V kn J 2 


( 64 ) 


+ L,.A. + iW(e«l-l) + i^(e“S 
V VW\ VW\ 

+ fw 2 ,/i 2 + 2 (e“" - 1)1 , 


1 ) 


where (•, •)2 denotes the inner product over Calculating the Gateaux derivative of 

I, denoted by dl = dl{w;w), we get 
2 


|2 
^* 112 


2 k 


11 


i=l 


\K\ 


+ (65) 


2 k 


11 


, , {k2l^A^\w^ + + - 1 + -h2 

K. I y fCii 2 

Lemma IV. 1 For the undetermined parameter fj, > 0 sufficiently large, there exists positive 
constants Ci, C 2 , and C 3 such that 


dl(w;w) - ^ ||Vwj ||2 > Cl 


2 = 1 


1 + lv^w^i 


■dx 


( 66 ) 


+ ^2 


^(k 2 i\Al^Wi + 11 ^ 1 ^ 2 )^ 

1 + T^|/C21 + \K\W2\ 

Kil 


■dx — C-i 


Proof. We note that the inner product terms of dl are of the general form, 


a ( V, e“°+^ -1 + ^1 , 

2 


for some positive constant a. Dehne 


M(n) = a — 1 + 


and let n = n’*' — n where = max{0,n} and v = max{0, —n}. So M{v) = M{v~^) + 
M(-n-). 

From the elementary inequalities, e* > a; + 1 and xy > — {x^ + y'^) for all x,y eM., we 
get 


M(n+) > ^ / {v^ydx-^ I {uo + gYdx>^ I {v^fdx-f3, 


a 


a 
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where /3 is a positive constant, since uq and g belong in For any a; > 0 we note that 


> 


X 


1 + X 


, therefore, 


M(n+) > 


a 




l + v 


-dx — (3. 


(67) 


Let ns now consider the eqnation 


M{—v ) = q;|^ —n , —1 + gj — g — ^ 

X 

From the elementary ineqnality 1 — e~^ > -for any x > 0, it follows that 

1 + X 

av~(l — g — ^ (1 — g) + a—^ -(1 — g — e“°). 


For p > 0 large enongh, we may obtain 1 — g > -. In addition, since 1 — e“° and g are in 
we have 


1 + f 


-{1 — g — e^°)dx 




< 


1 + X' 


-dx + / (1 — (/— e“°)^dx 


From the absolnte valne we can conclnde that 




-dx — (3, 


—-{1 - g - e'^°)dx >-\ [ '''' ’ 

1 + X 4 /to2 1 + X 


for some positive constant /9, since 1 — — e“° belongs in Thus, we have 

4 782 1 + X 

Consequently, fj671) and fl68l) gives us 

f 

M{v) = M{v*) + M(-ir) > Cl / , , , - Ci, 


( 68 ) 


(69) 


1 + |x| 

for some positive constants Ci and C 2 . Applying fl6^ to each inner product term of dl gives 
us the desired inequality fl66D . □ 

Absorbing constants into Ci and C 2 , we may rewrite (l66|) in the cleaner form 

2 . _ r 

1 ^YO) 


dl(w;w) - llVwill2 > Cl f 

i=i 


-dx 


+ C 2 


R2 (1 + I'iXll)^ 

{k2i\/\K\wi + |-A|tX2)^ 

(1 + \k2i\/\^\wi + |iC|x;2|) 


-dx — C 3 . 
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We reuse the positive constants Co,..., Cs and keep in mind that the importance of these 
constants is that they remain positive, of hnite value, and independent of wi and W 2 - 
Dehne the transformation T{wi,W 2 ) = {ui,U 2 ) such that 

ui = wi and U 2 = k 2 i^/\^\wi + \K\w 2 . (71) 


Applying this transformation to d70|) . we get 


dl{w] w) 


“ Je 


u; 


rdx + C 2 


Uo 


^dx — C 3 . 


(1 + ki)^ U{i + \u2\y 

Since T is an invertible transformation, there is a positive constant Co such that 


^l|VTO,||^>C„^||Vt„| 


2=1 


2=1 


Thus, we have the lower bound 

2 

dl{w-w) > CoY,\\^Ui\\l + Ci 
2=1 

which can be rewritten as 

2 


ui 


(l+|Ui|)2 


dx + C 2 


ul 


(1 + | u 2|)2 


dx — C 3 , 


dliw;w) > Co"r\\Vui\\l + (1 + \\Vu,\\l + [ 

i=i i=i V 




(1 + KI)" 

Using the following standard interpolation inequality over 

[ \v\‘^dx < 2 /" \v\‘^dx [ \Vv\‘^dx, v e 


dx I — C 2 . 


we may write 


vrdx < 4 


(l + |u|)^ 


dx / \v\‘^dx[l+ / \Vv\‘^dx 


Recall the generalized arithmetic and geometric inequality, n < < ^ AjOj, where the a^s 

2=1 2=1 
n 

and Ai’s are non-negative real numbers such that ^ A* = 1 . As a special case, we get 

2=1 


vl'^dx <C ( 1 -I- 


1 4 


(1 + |U|)2 

2 


dx 


1 4 


\Vv\‘^dx 


Ivrdx 
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It follows that, 


vl'^dx < cl 1 + 



n 

2 

(1 + 

|n|)2 


dx+ \Vv\‘^dx 


(72) 


for some positive constant C. By fl72|) . we obtain 


dl{w; w) > Co E l|V“j|li2(K2) + ||'Ui||L2(R2) - Cl. 


(73) 


i=l 


Using the relationship between Ui and Wi, we get the coercive lower bonnd 


dl{w]w) > Uo 1^11^111 1,2 + Ik 2 ||i, 2 ^ - Cl, (74) 

where Cq and Ci are positive constants independent of Wi,W 2 G lU^’^(M^). From (174)) . we 
conclude that, for any 5 > 0 there is an 0 < i? < oo so that 


inf dI(w,w)>S, (75) 

lhl|l,2=^? 

which gives the existence of an interior critical point of I in some open ball in H^. Therefore, 
the existence part of Theorem III. 21 is established. The uniqueness follows from the strict 
convexity of I. 


V. ASYMPTOTIC ANALYSIS AND QUANTIZED MAGNETIC FLUX 


In this section, we prove the exponential decay estimates of Theorem III. 21 We also 
establish the quantized magnetic flux integral formulas of Theorem III. 31 over the full plane. 


Our analysis follows similar ideas as in 


14 


,0,£i. 


Lemma V.l Let w = ( 101 , 102 ) in he the solution of the variational equations (16lap - 
(I61bp . Then Wi(x) ^ 0 as |a;| 00 , fori = 1,2. More explicitly, Wi satisfies the uniform 

decay 


lim sup |tai(a;)| = 0. 

|a;|=_R 


Proof. We will need the following embedding inequality 


38| for p > 2 


(76) 



P-2 

2p 


1 , 2 - 


(77) 
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Let us first justify that e-^ — 1 E for any / in By expanding {e-^ — 1)^ and 

using the Taylor series for e-^ we obtain 


J 


+5: 


k=3 


k\ 




(78) 


Integrating fl78|) over and using the embedding inequality fl77|l we get 

2*^-2 




< 


k=3 


k\ 


fdx 


(79) 


,^ 2^-2 


2 

2 ’ ^ k\ 

k=3 


k 


- 1 


fc -2 

2 


k 

1 , 2 ’ 


which is a convergent series. Therefore, e-f — 1 is in for any / in 

Let Wi, in be the solutions to (I61al) - (]61b[l . then they are also in for 

p > 2 by (ITT)) . In particular, the Wi are in To establish the desire uniform decay, 

it is sufficient to show that Wi belongs to ly^’^(M^). To this end, we rewrite the system 
fl61al) - fl61bD in the form 


Atci = 


2 K 


11 


(e“c. _ neV*!"! + - 1 ) 

VM 


\/W\ 


2 K 


(80) 


12 / u 


vw 

2Ku ^(\K\w2 +K2 i^\w,) _ 11 , 

vwr 

Aw 2 = 2{e< - 

+ - 1 ) + h 2 , 


(81) 


and show that the right hand sides of flSUD - flSTD are in L^(M^). 

We note that hi and /i 2 are in L^(M^). The terms, — 1 and — \ 

are in L^(M^) by fl79|) . Moreover, the term (e“o — l)eV^"'i is in L^(M^), by 

<m[ {e< - lydx < oo, 


where M = sup ^ ^ ^ L°°(R^) and e“o - 1 G Similarly, we 

a:eR 2 

can conclude that {e< - jg in L^(M^). Thus, the uniform decay of tci 

and W 2 is established. □ 
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Lemma V.2 Letw = {wi,W 2 ) be as stated in Lemma\V^ ThenVwi andVw 2 , also satisfy 
the uniform decay |Vwi|, |Vw 2 | —)■ 0 as |x| —)■ oo. 


Proof. We note that Wi belongs in for all p > 2 and z = 1, 2. We simply need to 

extend Wi to belong to for all p > 2. To achieve this, it is sufficient to show that 

the Awi is in for p > 2. This follows by proving that each term on the right hand side 

of (180])-(IHID is in Lp{R‘^) for p > 2. It is enough to show that ^(^e^\Kh^+k 2 iy/\K\wi)/ku _ 
and - 1) are in 

In general, if / is in for all p >2, we can conclude that 


|V(e^-l)pdx = j e^^Vfl^dx 

<m[ |V/pdx < cx). 


where M = sup < oo since / is in L°°{Rf) and V/ is in By fITTD . 


ceR2 


g(|J^k 2 +fc 2 i^«>i)/fcii_l are in Lp{R^) for allp > 2. The term (e“o-l)e^^"'l 

is in LPfR"^) for 2 < p < oo since 

< M J (e“o — l^dx < oo, 

where M = sup ^ qq since tCi G L°°(R‘^) and e“o — 1 G for p > 1. 

Moreover, (e“o — l)eV^"’i is in since |e“o — 1| < 1 and wi is in L°°{R^). Hence, 

{e< - is in Lp{R‘^) for p > 2. 

Similarly, we can show that (e“o —is in for p > 2. Therefore, 

the uniform decay of Vwi and Vta 2 follows. □ 

From {wi,W 2 ), defined in Lemma TV.ll we get (hi, 1 ) 2 ) by the transformation (IT0|) . and 
hence a solution pair (mi,M 2 ) is obtained as a solution of (I26ap - (l26b|l on the full plane M^, 
satisfying ([29]). To complete Theorem III. 21 we just need to establish the exponential decay 
estimates (jHOll and (l3lil . 


Lemma V.3 For the pair {ui,U 2 ) stated above, there holds the exponential decay estimate 


(mi + In 2)2 + {u 2 + In 2)2 < 


(83) 
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for |a;| sufficiently large, e G (0,1) arbitrary, > 0 a constant depending on e, and Aq = 
4min|2, 


Proof. The eigenvalues of the matrix K are Ai = 2 and A 2 = y. Both eigenvalues are 
greater than zero and since p q, we have that Ai 7 ^ A 2 . Hence, there is an orthogonal 2x2 
matrix O that diagonalizes K, i.e., 


O-^KO = 


Ai 0 
0 Ao 


:= A. 


Take Rq > max{|pi|,..., |gi|,..., IgAfjl}- Outside the disk of radius Rq centered at the 

origin, = {x G M^||x| < the system (I26ap - (l26bll considered over the full plane, 
becomes 

Aui = 4 {kiie^^ + ki 2 e^^ — 1) where i = 1,2. (84) 


Rewrite (|84|l in the form. 


Aui = AkiiUi + Aki2U2 + 2/cii(2e“^ — 2 mi — 1) + 2/ci2(2e“^ — 2u2 — 1). (85) 


We would like the terms 2e“‘ — 2 m j — 1 to converge to 0 as |x| 00 . Thus we dehne 

Mj = Mj + ln(2). It then follows that each Xj —)■ 0 as |x| —)■ cxo in the sense of Lemma fV.ll 
Then fl85l) becomes 


Avi = 2kiiVi + 2ki2V2 + 2Kii{e^^ - xi - 1) + 2ki2{e''^ -V 2 -I). (86) 

Dehne the new variables Ui and U 2 such that 



Using the variables Ui and U 2 , we express in the form 

/ Af/, \ ^ / 2X,U, - ^(U, + U2)-l\ 

[AU 2 J \ 2 X 21/2 J - ^([R - 1 / 2 )-ij 

By using the Taylor expansion of e*, we may write 


AR, 


2 XiUi + anfUi, 112)1/1 + ai 2 {Ui, 1/2)11/2, 
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where ajk{Ui, U 2 ) —?• 0 as |x| —)■ 00 (j, k = 1,2). In two dimensions, we recall the inequality 

Ay = 2/A/ + 2(|)%2(|)%2/A/^ 

The following inequality then follows, 

+ Ul) > \oiUl + Ul) - aiU^, U2){Ul + Ul), 


where the function a{Ui, U 2 ) —)• 0 as |a:| 00 . Consequently, for any e G (0,1) we can find 

an i? > i?o large enough, so that 


A(17/ + Ul) > (1 - 1) HUl + Ul), |a;| > R. 

(87) 

Introduce the following comparison function. 



^(x) = |a:| > 0 , 

C, a G M, C, cr > 0. 

( 88 ) 

Then, 

X 


Ae = - 

(89) 


Subtracting fl87|) and flH^ . for | x\ > R, we get 

A(t/? + C/J - e) = A(t;f + ul)- A? 



> (1 - 0 + UD - 


We have the freedom to select cr^ = (l — |) Aq, which gives us 

A{Ul + Ul-i)>a'^{Ul + Ul-i) for \x\ > R. (90) 

Now we can select C in fl 88 l) large enough so that Uf + U 2 — ^ < 0 for |x| = R. Let’s denote 
C by to point out its dependence on e. We would like to extend this inequality so that it 
holds for all |x| > R. Since f/f + f/| —)■ 0 as |a;| —)■ 00 and applying the maximum principle, 
we can conclude that t/f + Cf — ^ <0 for |x| > R. For any e G (0,1) we have the following 
useful inequality y/f — e/2 > 1 — e. Hence, for |a;| > R we get 

Uf + ul < C',e"Vd-^/2)^oh| < 

By the orthogonality of O, the equation Ul + Ul = v\ R follows. Therefore, for |a:| > R, 
we have the desire inequality 

vl + vl = (mi + In 2f + {u 2 + In 2)2 < . □ 
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Lemma V.4 ui and U 2 , from Lemma W.^ also satisfy the inequality 


+ iVwsr < 

where e, C^, and Aq are as defined in Lemma fK At 

Proof. Take R as defined in Lemma [V.3[ Differentiating equations fl 86 |) . for |a:| > R, we 
get 

AdjVi = 2kiidjVi + 2ki2djV2 + 2kii{e"^ - l)djVi + 2ki2{e'’^ - l)djV2., (91) 

for i,j = 1,2 and dj = 

Define U = {Ui,U 2 Y = ( 91 ^ 192 ^^ 2 )"^ and E{x) = diag{e“i*^*\ where diag{a, 6 } is a 

2x2 diagonal matrix with diagonal entries a and b, respectively. Then the system (l?I|) may 
be rewritten in the form 

AU = 2KU + 2K{E{x) - RW, 

where R is the 2 x 2 identity matrix. Consequently, we can establish the following inequality 
for |a;| > R, 

A\U\^>Xo\U\^-b{U,,U2)\U\\ 

where b{Ui,U 2 ) —)■ 0 as |a;| —)■ 00 . Hence, as in Lemma IV.31 for |a;| > R, we arrived at the 
inequality 

+ f/| = {djViY + {djV2Y < 

where are positive constant depending on e. Therefore, we can take C^ = 2 maxjCg^i, C^^ 2 } 
and obtain the desire inequality 

|VmiP + |Vu 2 p < for |x| > R. □ 


As a result of the exponential decay estimates, we get the quantized magnetic flux inte¬ 
grals. A direct calculation shows that the integrals of the functions 5 'o(x) and g'Y'^) over 
are independent of the parameter p. More explicitly. 



47 rA^i and 



= 47rA'"2. 


The divergence theorem in two dimensions gives 


Av^dx 


Av2dx = 0 . 
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By integrating equations fl60al) - fl60bD over the full plane, we get 



+ 


l)dx 


-vriV,. 


Similarly to the doubly periodic case, we obtain the quantized magnetic flux integral 
formulas of Theorem III. 31 



—2npNi 


and 



—2'KpN2. 
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